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Abstract 

Using a recently developed perturbation formalism based on curvature 
quantities, we complete our investigation of the linear stability of black holes 
and solitons with Yang-Mills hair and a negative cosmological constant. We 
show that those solutions which have no linear instabilities under odd- and 
even-parity spherically symmetric perturbations remain stable under even- 
parity, linear, non-spherically symmetric perturbations. Together with the 
result from a previous work, we have therefore established the existence of 
stable hairy black holes and solitons with anti-de Sitter asymptotic. 



1 Introduction 



Soliton and hairy black hole solutions of general relativity coupled to non- 
Abelian gauge fields have been the subject of intensive investigation for over 
ten years (see [||] for a comprehensive review of the subject). The field was 
sparked by the discovery of solitonic [Q] and coloured black hole |^] solutions to 
su(2) Einstein- Yang-Mills (EYM) theory. Since then, many examples of both 
solitons and black holes in various theories involving non-Abelian gauge fields 
have been discovered, both with and without a cosmological constant. 

Many of these examples have been shown to be classically unstable, making 
the search for stable solitons and hairy black holes vital in this area. An im- 
portant (and surprising) discovery was that there are solitons and black holes in 
su(2) EYM with a negative cosmological constant which are stable with respect 
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to linear, spherically symmetric, perturbations |[. This result was surprising 
because all such solutions in EYM with either a positive or zero cosmological 
constant must be unstable |^. For those solutions which are stable with 
respect to spherically symmetric perturbations, it remains to be seen whether 
the stability is altered if non-spherically symmetric perturbations are considered. 
The perturbations fall into two sectors, depending on their behaviour under par- 
ity transformations. In Q we analyzed the odd-parity sector of perturbations 
and showed that those solutions which are stable under spherically symmetric 
perturbations remain stable under odd-parity, linear perturbations. The proof of 
stability will therefore be complete if we can show that there are no instabilities 
in the even-parity sector. This is the subject of the present article. 

In common with we shall use a recently developed perturbation formalism 
based on curvature quantities. The main advantage of this formalism IC] is 
that on a static and purely magnetic background it yields a wave equation for the 
linearized extrinsic curvature and electric YM field where the spatial part of the 
wave operator is (formally) self-adjoint. This offers the possibility of studying the 
linear stability of non-rotating black holes and solitons by analytical means. For 
example, we can use the nodal theorem |11] which in fact will turn out to be an 
essential tool in our stability proof. Even with these powerful tools at hand, the 
even-parity sector is considerably less amenable to analysis than the odd-parity 
sector. As we have shown in our previous work the stability in the odd- 
parity sector is topological in the sense that we do not need to know the exact 
details of the background solutions in order to show that the system is stable. 
In the even-parity sector, we will see that the detailed structure of the metric 
is needed even in some simple cases. Nonetheless, we will present an analytic 
proof that those solitons and black holes having no instabilities under spherically 
symmetric perturbations remain stable under even-parity linear perturbations. 
In short, all soliton and black hole solutions when the cosmological constant is 
sufficiently large and negative are linearly stable. 

This work is organized as follows. In section ^ we remember some impor- 
tant results about the hairy black hole solutions with a negative cosmological 
constant, recently found in [^, and the corresponding solitonic solutions |^. 
In section ^, we briefly review the curvature-based formalism of perturbation 
theory for a static background. The harmonic decomposition is performed in 
section ^ Special cases, such as the stability of the Schwarzschild-anti-de Sitter 
and the Reissner-Nordstrom-anti-de Sitter metric are discussed in section ^. In 
section |6|, we discuss the general case and prove that the system is stable when 
|A| is large enough. Given that the proof is rather involved, we first present 
an outline of our argument in section 6.1 before working through each step in 
detail. Technical details such as the discussion of the perturbation potential in 
the Reissner-Nordstrom-anti-de Sitter case and the factorization of the spatial 
perturbation operator in the odd-parity sector are discussed in appendices ^ 
and respectively. 

The metric signature is (— , -|-, -|-, -|-) throughout, and we use the standard 
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notations 2a;(afe) = ujab + ^^ba and 2u;[ab] = ^ab — ^ba for symmetrizing and anti- 
symmetrizing, respectively. Throughout the paper, greek letters denote space- 
time indices taking values in (0, 1,2,3), while roman letters will denote spatial 
indices taking values (1,2,3). 

2 Solitons and hairy black holes with a negative cos- 
mological constant 

In this section we discuss very briefly those details of the spherically symmetric 
black hole and soliton solutions of su(2) Einstein- Yang-Mills theory with a neg- 
ative cosmological constant that we require later in this paper. Further details 
can be found in |^] and the original papers [Q, ^. 
The equilibrium metric is spherically symmetric 

ds'^ = -N{r)S^ir) dt'^ + N-\r) dr"^ + {dO^ + sin^ 6 d(t>'^) , 

and the gauge field potential has the spherically symmetric form 

^ = (1 - w{r)) [-T^de + Te sine dcj)] . 

Here the su(2) generators Tj.^e,(/) are given in terms of the usual Pauli matrices 
cTj by Tr = ■ g^/2i, etc. Writing N{r) = 1 — 2m{r)/r — Ar'^/3, where A is the 
(negative) cosmological constant, the field equations take the form: 

S 'p 


where G is Newton's constant, and we have set the gauge coupling constant equal 
to ^/Att for convenience. Here, and in the rest of the paper, we use a comma in 
a subscript to denote partial differentiation, i.e. w^r = drW. 

For solutions which approach anti-de Sitter (adS) space at infinity, the asymp- 
totic behaviour of the field functions is: 

m(r) = M + — + 0(r"^), 
r 

w{r) = Woo-\ \-U{r ), 

r 

S{r) = l + 0(r-^). 

Due to these boundary conditions, in general the solutions will be globally mag- 
netically charged. 
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For black hole solutions having a regular event horizon at r = r^, all the field 
variables have regular Taylor expansions near the event horizon, for example 



w{r) = w{rh) + w^rirh){r - Th) + 0{r - r^f . 

However, there are just two independent parameters, S{rh) and w{rh) since 
= at the event horizon, which gives 

■m{rh) = — - — . 

In order for the event horizon to be regular, we shall also require that N^ri'i'h) > 0, 
which implies that 

f.^i-a.i-gM:^>o. 

From (|), one has 

'"''^''^ = ■ 

There are also globally regular (solitonic) solutions, for which the behaviour near 
the origin is: 

m{r) = 2Gb'^r^ + 0{r^), 
w{r) = 1 - br'^ + 0{r^), 
S{r) = 5(0) [1 + + 0{r^)] . 

Here the independent parameters are b and S{0). 

The simplest solutions to the field equations (|l|) are the Schwarzschild-adS 
solution, 

w = ±1, 5=1, m = const, 
and the Reissner-Nordstrom-adS (RN-adS) solution 

G 

w = 0, 5 = 1, m = const. . 

2r 

In both cases, the YM field is effectively Abelian. It is however interesting to 
study the stability of these solutions with respect to non- Abelian perturbations 
(see section |5|) . 

The solutions of greatest interest in this article are those effectively non- 
Abelian solutions, for which the gauge function w has no zeros, since these 
solutions were shown in Q to be linearly stable to both even- and odd-parity 
spherical perturbations. In Q it is proved that for any value of the gauge field 
at the event horizon, w{r}i) ^ 0, for all sufficiently large |A| there is a black 
hole solution in which w has no zeros. Similar behaviour is found numerically 
for the solitonic solutions |5|. For spherical perturbations of these equilibrium 
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configurations, it is proved analytically that all the solutions in which w has 
no zeros are stable in the odd-parity sector Q. The even-parity sector is more 
complicated, but stability can be proven for sufficiently large |A|. 

In addition to the boundary conditions on the field variables at the origin (or 
event horizon, as applicable) and infinity, we shall in section ^ make extensive 
use of the fact that the equilibrium solutions (both solitonic and black hole) 
are analytic in both r and the parameter A. This is proved in |^ for black 
hole solutions, a proof which is readily extended to cover the solitonic case (see, 
for example ||l2| for the asymptotically flat situation). We shall also apply the 
result, proven in H, that 
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w^r{r) ^ o{\A\ 2) as |A| ^ 00. 

Dyonic solutions with a non-vanishing electric field also exist in this model 
with A < ||]. However, in common with our analysis of the odd-parity per- 
turbations [g], our work here applies only to the purely magnetic equilibrium 
solutions. In addition, recently hairy black holes with non-spherical event hori- 
zon topology have been found in this model |]13[. Here we consider only black 
holes whose event horizon has spherical topology, although we conjecture that 
all black holes (whatever the topology of their event horizon) will be stable if 
|A| is sufficiently large. 



3 The pulsation equations 

Generalizing previous results ^ , we have shown in |8| that linear fluctuations 
on a static and purely magnetic background are governed by a symmetric wave 
equation for the linearized extrinsic curvature and electric YM field. Performing 
an ADM decomposition of the metric and the gauge potential, 

g = -a^dt^ +gij{dx' + /3'dt){dx^ + P^dt), 
A = -^adt + Ai{dx' + P'dt). 

(where $ and Ai are both Lie algebra valued), these quantities are given by 

Lij = ^6 {dtQij - Lfsgij) 

and 

Si = -6{dtA + Di{a<^)). 

(We use the same notation as in In particular, a bar refers to the 3- metric 
gij or to the magnetic potential Ai .) On the background, static coordinates 
are chosen such that the shift /3* and the time derivative of the 3-metric, dtgij 
vanish. Similarly, the gauge is chosen such that the electric potential <I> and the 
time-derivative of the magnetic potential, dtAi vanish. As a consequence, Lij 
and £i are vector- invariant, that is, invariant with respect to both infinitesimal 
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coordinate transformations within the ADM slices and infinitesimal gauge trans- 
formations of the gauge potential. More precisely, with respect to a coordinate 
transformation 5x'^ 6x^ + X^^, we have 

Lij ^ Lij + V^i{a^V^)f), (2) 
Si ^ £i + a^F^jV^f, (3) 

where / = X*. The fact that the spatial components, X*, do not appear in (^) 
and (|3|) justifies the name "vector-invariant" . 

The pulsation equations can be obtained from the following energy func- 
tional: 

E = Egrav + Eym + Eint , 

where 

If/ Q^j^i 

Egrav ~ 2 ( Q,2 ^ij-^ki + G^'' C^'^ Lij){V gLki) + 2L'^^ R ^Ljk — R^uj L 
- 2kV^U, + 4GTr U'^F^P^Lki + ^FmF^'L,,lA\ a^dx\ 



Eym = G I Tr{ ^S'Si + 2D^'£^^ ■ D^iSj] + 



( ^ 

a 



Eint = -2G I Tr <j VWk ( aFi'£, ] - aF^'rVkU, + -V'£kD, ( a'F,' 



a 



- 2F''\Vka)L£i + aFij£'V^ ^dx^. 
Here, A is defined by 

A = L, 4 

a 

which is the difference between a vector-invariant combination of the time- 
derivative of the lapse and the trace, L = g^^ Lij , of the extrinsic curvature. 
By absorbing L into A, one makes sure that the resulting dynamical equations 
are hyperbolic. In fact, it can be checked that (^) corresponds to a linearized ver- 
sion of a densitized lapse, which is widely used in recent hyperbolic formulations 
of Einstein's equations (see Ref. |jl^ for example). 

In terms of the variable A, the harmonic gauge reads A = 0. In this case, 
one has a set of equations with symmetric potential, but the kinetic energy is 
not positive This is the reason why the De Witt metric 
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grav 



appears in E, 

As we have discussed in [10], one can also adopt the maximal gauge, L = 
0, in which case the equations are both hyperbolic and symmetric, but where 
perturbations of the lapse still appear in A. In this case, one recovers the energy 
functional given in Appendix A of Ref. Q. At this point we recall that for 
perturbations with odd parity, L = and A = anyway since they represent 
scalar quantities. In the even-parity sector, it will turn out to be useful not to 
choose a particular gauge a priori but to derive the pulsation equations in a 
general gauge and to choose an appropriate gauge later. 

The dynamical equations resulting from the variation of E are subject to the 
linearized momentum and Gauss constraint equations, 

= Ci = aG/^'Vj - 2GTr (F/,- 

= g = a& . 

Additional constraints involving also perturbations of the metric and the gauge 
potential themselves are the Hamiltonian constraint and all evolution equations, 
which we had differentiated in time in order to construct the wave operator. 

Finally, we will make use of the following fact later in our stability 
analysis: In any gauge, the terms involving A in the energy functional do not 
contribute provided that Lij and £i satisfy the momentum constraint Cj = 0. In- 
deed, using partial integration (and homogeneous Dirichlet boundary conditions 
for the perturbations), the terms involving A are 

/ {^-]^^'i^^Lij^ka^^iA - 2GaTr \^Fij8'V^ A^^ ^dx^ 

= j a [aVk [g''''^] - 2GTr {F^'Sk}^ ViA^gdx^ 

which vanishes if C' = 0. 



4 Even- parity fluctuations 

Here, we specialize the pulsation equations given in the last section to a spher- 
ically symmetric background with gauge group su(2). In this case, it is conve- 
nient to expand the linearized extrinsic curvature and electric YM field in terms 
of spherical tensor harmonics since then, perturbations belonging to different 
choices of the angular momentum numbers i and m decouple. Furthermore, the 
tensor harmonics can be divided into parities. Odd-parity perturbations were 
discussed in |^]. In this work, we consider the even-parity sector. 
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4.1 The vacuum case 



We start by computing Egrav ■ For technical reasons, it is convenient to parametrize 
the background 3-metric according to 

g = dx^ + r{xfdn'^, 

where x is a radial coordinate, and where r and the lapse a depend on x only. 
In terms of the functions and S defined in section |2|, we have a'^ = NS"^ and 
N = {d^rf. 

Einstein's background vacuum equations are obtained from 

- Roo = ^ [w H —5 - A, 



-2 = Rxx = ^ w'^ - ^ + A, 



2(l-r'^-rT" -rr'-\= R\ =^i^^^!^ + 2A. (5) 



Here and in the following, capital indices refer to coordinates on the 2-sphere, 
and a prime denotes differentiation with respect to x. 
In the even-parity sector, we expand Lij according to 



T ~ (2) 

LxB = , 

Lab= (-^P + ^ij e% + g e%, 



r 



= rgABY, e% = r I^Va^bY + ^{1 + l)gABY j , 

form a basis of even-parity tensor harmonics, which are orthogonal with respect 
to the inner product induced by g. {Y = Y^"^ denote the standard spherical 
harmonics. The indices im are suppressed in what follows.) Also, we have 
chosen the parametrization such that t and p correspond to the trace and the 
radial trace-less part, respectively, of Lij. After the rescaling 



where 



with = £(^£ + 1)^ A = /x^ - 2, the expansion is normahzed such that the De 
Witt metric becomes 

The signs reflect the signature of the De Witt metric. 

Inserting the expressions above, using the background quantities 



^ AxB = -rr"gAB, R'hAB = 2(1 - n6^A9B]C, 
Rab = {l-r'^ -rr")gAB, 



Rxa 



r 



FxB = -w'i^TA , Fab = (w - l)Tr eab , 



and integrating over the spherical variables, the gravitational energy functional 
becomes, after some calculations. 
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grav 



{V , TV) + {dpV , TdpV) + {V , SgravV) - 2{V , r6(a) 



J grav / 



dp, 



where V = {t,p,q, g)^ ■ The matrix T is given by 





T -- 


= diag(- 


1,1,1,1) 




and the symmetric matrix Sgrav 


by 








( Stt 


Sti 





\ 


^grav ~ 


Sti 
^ 


Sii 


'S'22 




2VA7,p 
S'33 / 




^ 








with 
Sti 
Stt 



S22 
S33 



-Y(7«')' + -^7' 



w 



/2 



(W^ - 1) 



21 



14G 2 



2 13 , a' 5 4 ,5 , 2 

+ rr" H rr' \ a" ^a'^ + Ar^ 

3 a 6 a 3 



2r2 



^11 = 



7 



4G 2 



,, 31 .ol 4r2 

rr ' 

3 

iw' - l)2n 



/i^ + 6r'2 - 5rr" - — rr'— + --^(aa')' - 2Ar2 
3 a 3 



2r2 



7 



7 



/ 2 



^2 _ 2/2 _ ^^// _ ^^/^ _ 2Ar2 + 4Gu;'2 _ 2G 



,(^'-l)' 



a 
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Here, we have also introduced ^ = a/r and dp = adx ■ The inhomogeneous term 
b{a)grav ) where a = a(x) is the scalar amplitude parametrizing A, 



is given by b{a)grav = 

bt{a) 

61(a) 
62(a) 
63(a) 



A = aY, 
bt,bi,b2,b3){d), where 



,a 



2 1 



r(a^a')' + 2a^r'd' - n^— d 

r 



r{a d! 



1 1 
2 , , 

a r a H a 

2 r 



'Ifi"^ ar [—d 



a 



The resulting perturbation equations are 

(Ta/ - Td^ + Sgrav) V - Th{l 



Igrav 



0. 



(6) 



The linearized momentum constraint yields the following two equations: 







= aCT 



ar 



t 

ar 



„ r'^q 



2^/3/i 



It 



2 



■79 



-7P- 



r 



75 



Note that the spatial operator which acts on V is symmetric (with respect to 
the standard L? scalar product). However, compared to the odd-parity case, the 
following problems arise here: First, the kinetic energy is not positive. In order 
to fix this, one could multiply equation ^ from the left with T. However, one 
would then lose the symmetry of the spatial operator. Next, the (densitized) 
lapse still appears into the equations, and one has no evolution equation for d. 
On the other hand, the amplitudes are not fully coordinate-invariant but are 
subject to the gauge transformations (^). As a consequence, we have 



V + 6(0 



grav 



(7) 



where ^(x) parameterizes / according to / = ^y. Since the same operator- 
valued vector b(.)grav appears in both the coordinate transformation and the 
terms involving d in the dynamical equations, it is clear that any gauge-invariant 
combination of the equations for t, p, q and g annihilates the terms which depend 
on a. 
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Using the background equations (|5|), we can re-express rr", r'a' and a" in 
terms of = r'^ and matter terms: 







-N)- -r^ 


la 

rr — = 

a 






2"" 

a 


-(1 


( 

N) + 2G^ 



.12 



- 1)^ 



As a result, oiiG Ccin rewrite the coefficients of Sq^av 



Sti 
Stt 
Sii 

5*22 

S33 



4a^^ 272 
3 a 3 



1-N -Ar^ + 4Gw'^ - 3G 



X + 1 + N - Ar' - iGw'"^ + 3G 



A - 5 + 13iV + 5Ar2 - AGw''^ + ^ 



— Sti 



2 



(rO. 



a \r^ / ,pp 
r 



A + 4G 



- 1)2 



(8) 



4.2 The pure YM case 

Next, we compute Eym- The background gauge potential has the components 

= 0, Ab = {1- w)e%Tc , 
where w depends on x only. The background YM equation reads 

I , w"^ -I 



— (au;') 



w- 



a 



The electric YM field is expanded into su(2)-valued spherical harmonic one-forms 
with even parity (these one- forms are explained in Appendix D of Ref. p!5[|): 













r 




— ^ B 



1 



CTrVAY + dYTA + e + ^I^ Yta 

where in terms of the Pauli matrices a = (fx*), r,. = ■ a/ {2i), ta = Cj^- g_l (2i). 
Here X2 and X3 are the su(2)-valued harmonics 
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After the rescaling 

~ b _c ~ d ^ e 

the amplitudes are normahzed such that 

2G j Tr {g^^£i£j) y/^dx^ = j (b"^ + + d^ + e^)dx. 
The YM energy functional becomes 

Eym = \j {{W,W) + {dpW , dpW) + {W , SymW)') dp, 

where W = (6, c, d, e)"^ and 

^ + ^ + / + ^ sym. sym. sym. ^ 



7^ 



7 

-:^(7'«^),P - V A + 2/ + ?;2 sym. sym. 



V2/x^ ^/2/i«; A - 2 + 3/ + ^ 



sym. 

2Aif -\/2Aw; ^u^ _ 



(9) 
with 

f^w^ + l, n^2VG^, ^;^2^/G^^^. (10) 

(Note that the u here differs from the u in our previous article by a factor of 
2\/G.) The linearized Gauss constraint gives 

= V2Gfiarg = i^-fdp (^^^ + -fwc - ^ - y^Tej X3 . 
4.3 The interaction term 

We finally compute the interaction energy. Using partial integration and the 
background equations D^[aFjk) = 0, it is convenient to rewrite Eint as 

Eint = -2G J |2L^^Tr (^F,''{D£)kj^ + 2vHi + 2^ VkL + v'^VkA^ a^dx\ 

where k = aV^{Lij/a) and = Th:{F^^£k)- 

Inserting the above expansions for Lij and £i and integrating over the spher- 
ical variables, we obtain 

Eint = J {{V, A^ntdpW) - {AltdpV , W) 

+ {V, SintW)-{bint{d),W))dp, 
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where 



Hnt 



U 

2 







4 



'int 





1 
V 

_ r fr 











V2 J 



V6 



+ 



^ 

■q U ^ — 



7 \ 



V2 



71 / 



(11) 



and 



a 



b'intia) = [^ua,-n^va,V2rua',0] . 



At this point, we remind the reader that although the hnearized electric 
YM field is invariant under infinitesimal 5u{2)-gauge transformations, it is not 
invariant under infinitesimal coordinate transformations, see @. As a result, it 
is easy to show that the YM amplitudes W transform according to 



W^W + bintiO, 



where f = . 



4.4 Summary 

Taking together the above results, the total energy functional becomes 
E = ]^U{U,TU) + {dpU , TdpU) + {U , AdpU) + {U , dpAU) 



+ {U, SU)-2{U,Tb{a)grav))dp, 
and the pulsation equations are 

{Tdt^ - Td^ + Adp + dpA + S- Tb{a)) U = 0, 
where U = {V, Wf, T = diag(-l, 1, 1, 1, 1, 1, 1, 1), 



(12) 



(13) 







, s 



grav iJint 



and 



b{d) = {bgrav{a),bint{a)) . 
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Infinitesimal coordinate transformations are given by 



U^U + biO, (14) 
and the constraint equations are 

= ^2dp{^)-»rd,^-^y^jq-^nd, (15) 

= -^Opi 5— 7* — 7p- VA7sr-n6 + 7t;c, (16) 

\ a ) 3 3 

= 79p (^^) +7'«^c- -^7^- y^7e. (17) 

Finally, as in our previous work, we assume that all perturbations vanish at the 
boundary points, so that all boundary terms vanish. 



5 Special cases 

In this section, we discuss the perturbation equations for some special cases. 
We first show that for £ = 0, one recovers the radial pulsation equation derived 
in Q. Then, we introduce fully gauge-invariant amplitudes and discuss the 
stability of the Schwarzschild-adS and the RN-adS solutions. These examples, 
which are much simpler than the cases where the gauge potential is effectively 
non-Abelian, indicate that a topological stability analysis, as performed in the 
odd-parity sector, is not possible in the even-parity sector. 



5.1 Radial perturbations 

For £ = 0, the only gravitational amplitudes are t and p. Since these are sub- 
ject to the one-parameter family of coordinate transformations ( |T^ ) and to the 



constraint equation (15), one expects that there are no dynamical gravitational 
modes. In fact, one can show that for £ = 0, the only physical gravitational 
perturbations correspond to the variation of the mass. This is also clear in view 
of Birkhoff 's theorem. For the YM field, we have only the amplitude d and the 
Gauss constraint is void. Therefore, we expect to have an unconstrained wave 
equation for the YM field when £ = 0. Looking at (0) we see that we can 
construct a gauge-invariant linear combination: 

; , \ ru , , 
(i = (i + ^— -t+p). 

For computation, it is convenient to choose a gauge in which t = p since in this 
gauge, d = d. Furthermore, in this gauge, the only constraint equation yields 
the simple relation 

ru 
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Using this, and taking the combination of the dynamical equations which corre- 
sponds to the definition of d, we obtain the following wave equation for d: 



{d,'-d^)d + 2 



ru 



u,p d 



1 + 



ru 



Finally, using ru/r^p = 



-1 + 2N + Kr^ + G 



d = 0. 



AG'j'^Nw^ and the background YM 



equation ru^p = 2^/G'^'^w{tiP' — 1) — ur^p , we have 



l))d 



+ 4G72 



2ww , 



w 



w , 



l-Ar^ -G 



- If 



d = 0, (18) 



which agrees with the equation found in [Q] . By replacing d by d, this equation 
acquires a gauge- invariant meaning. Using the estimate of section |^ for A — > — 00, 
we see that the dominant terms in the potential are 7^(3t(;^ — 1) which is positive 
for solutions with w{rh) > l/\/3 and large |A|. This shows the existence of 
solutions which are stable with respect to even-parity radial perturbations @. 

5.2 The gauge-invariant approach 

Motivated by the analysis above, one can try to introduce gauge-invariant am- 
plitudes for I >1. In order to do so, it turns out to be convenient to introduce 
the amplitude 

Using the transformation rules (|^), we find 

a2 

T T H ^, 

r 



where 



fr = — + rr — 
Z a 



A 3m „ 
N = - + — + Gw' 
1 r 



/2 



G(w^ 



I) 



2 

The advantage of the new amplitude r is that - like the amplitude g - it trans- 
forms with an additive term which is algebraic in ^. Provided that has no 
zeros, one can use r in order to construct gauge-invariants. For example, we can 
construct 



( = 9 



2 fr 



(19) 



which is gauge-invariant for £ > 2. In fact, we will show below that for vacuum 
gravity, <^ satisfies the Zerilli equation. 
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Using the momentum constraint equations ( [TsD and (|T6|), we find the follow- 
ing constraint equation 



7 



+ 



X r ri r'u , 1 r u ^ , ^ 

r- ' w g + ^^b-^ — '■^vc+-d. (20 



In the gauge where r = 0, this yields an algebraic relation between q, g, b, c and 
d. 

Next, we compute an evolution equation for r: Taking the combination of 
the evolution equations which corresponds to the definition of r, we arrive at 
the equation 

2-1)2- 



□ r + 7^ 






6m 








r 






+ 7 







a 



,,'2 



T + ^ijr 

J \ 7 
(u;2-l)2- 



-yuvw 



J V2fx 



[w 



If 



+ 



1 



r'udpb + 2r'unb + r'— ub + ;U^7n 6 + 2ar"ub 



1 

+- i^ne fra, 

IjjL r r 



I rr Y \ o " 
jwr u — V \ — 2ar jv 

1 \ r 



a 



c + I u 



(21) 



where we have defined 

We have also used the momentum constraint equation (|T6|) in order to eliminate 
first derivatives of q. It seems clear from equation (^) that a formulation of the 
evolution equations in terms of gauge-invariant quantities would be quite messy. 
Worse than this, we have no reasons to expect that the resulting equations have a 
symmetric potential. Nevertheless, in some special cases where w is constant, the 
expressions simplify enough and one can find a fully gauge-invariant description. 
This will be the subject of the next two subsections. 



5.3 Stability of the Schwarzschild-adS solution 

If {wl = 1 (i.e. the YM field is a pure gauge), we see that the interaction terms 
Aint and Smt vanish, and the gravitational and YM perturbations decouple 
from each other. In order to investigate the gravitational sector, we adopt the 
gauge-invariant formalism described in the previous subsection. First, we note 
that for £ = 1, the constraint equation (pO|) yields, in the gauge r = 0, 



q = 0. 
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Then, the evolution equation for r, (pTj), gives a = 0. This shows that for 
the gravitational sector is empty. 

For £ > 2, we introduce the gauge-invariant amplitude 



1, 



C = 9 



H^X 1 



2 fr ' 

where = A/2-1- 3m/r. We derive the equation for ( in the gauge r = 0. 
Taking the combination of the evolution equation for r, (21), and the equation 
for g corresponding to the definition of C, the terms involving a cancel, and we 
obtain 



^ \/A 9 ,6m 
2 r/ 



S33-7'^(l-2iV-r2A) 



Next, the constraint equation ([20|) yields 



Using this, and r 



N 



g = 0. 



q = -y/X^g. 

It 

1 — 2m /r — Ar'^/S, we eventually get the Zerilli 



equation ||T^ (with the presence of a cosmological constant) 

{n + Vz)g = Q, 



(22) 



where 



Vz = N 



A^H[(A + 2)r + 6m] + 12m^(3Ar + Qm - 2AH 



(Ar -|- 6m)^r^ 

If we replace g by this equation becomes manifestly gauge- invariant. Since 
Vz is positive for all r G (r^, cxd), there are no gravitational instabilities. At this 
point, let us also mention that equation ( [2^ ) and its odd-parity counterpart have 
already been derived and used in order to compute the quasi-normal modes of 



Schwarzschild-adS black holes (see Ref. |17] and references therein). 

The YM perturbations describe a linearized YM field on a Schwarzschild 
black hole. We will show later (see section p.2| ) that those fields cannot possess 
exponentially growing modes either. 



5.4 On the stability of the RN-adS solution 

When w = 0, we see that the equations decouple into two sets, one set comprising 
the amplitudes {t,p,q, g,c) (referred to as Abelian amplitudes in the following) 
and the other set comprising (6, d, e) (referred to as the non- Abelian amplitudes). 

We start with the discussion of the Abelian amplitudes: The relevant gauge- 
invariants are 



c 



g - 

c + 



fi^X 1 
2 // 

/i V 
V2Tr^- 
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(For i = 1, is void.) Taking the combination of the evolution equations cor- 
responding to these ampUtudes and ehminating q using the constraint equation 
(pO|), i.e. using 



r' 

Ir L 



^9 



vc 



c 



0, 



we get the following system of equations: 

where the potential Vm is symmetric and has the particular form 

U + W 



(23) 



M 



N 



-3M VaGX 
y/iGX 3M 



Using N = 1- 2M/r + G/r'^ - Kr'^/2> and 2/^ = A + 6M/r - AG/r"^, we can 
write the functions U and W as (the fact that > for r > r/^ is shown in 
Appendix ^) 



U 



1 



r2(2^)2 
+ 4 



+ 2 + 



9M 4G\ 



J.2 j 



2 ^ . 3M 9M2 + 2G 
-GA + + 5 



4GMA 9M2 
+ - 




16GM 6G2 
^ + —r 



32AfG2 



8G3 



W 



A^ + 4A + 



AM ^ G^^ _^ 2Ar2 ^6M BGy 

r 



+ 



3 \ r y 



For A = 0, these equations reduce to the ones obtained by Moncrief |18] de- 
scribing even-parity perturbations of a RN black hole. The fact that Vm can be 
written as the sum of a function and a function times a constant matrix allows 
the system to be diagonalized. The absence of exponentially growing modes in 
this system is discussed in Appendix |A|. 

For ^ = 0, the non-Abelian modes are described by equation (18) where one 
sets w = 0: 

(□ -7^) d = 0. 

In |§, we have shown that this equation possesses exponentially growing modes. 
Therefore, the RN-adS solution is unstable with respect to non-Abelian odd- 
and even-parity radial perturbations. 
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For £ > 1, there are no instabihties. In order to see this, we consider non- 
Abehan perturbations, which are described by the system 







□ 6 + 



^ + 7'(A + 1) 

7 



(For . 



6 + \/2//7,p(i + \/2A7,pe, 

= □d + \/2/i7,p6 + 72(A + l)d, 
= □e + \/2A7,p6 + 72(A + l)e. 

1, the equation for e is not present.) We also have the Gauss constraint. 



= ldp 



V2 



7 d ■ 



7 e. 



Using this constraint, we can eliminate both d and e in the evolution equation 
for b. Defining B = b/'j, one obtains the equation 

nB + ^\X + l)B = 0, 

which has no unstable modes. From the Gauss constraint, it follows that ^d + 
e cannot grow exponentially in time either. Finally, we can find an equation 
for F = d/n- e/\/A: 

□ F + 72(A + 1)F = 0, 
which is identical to the equation for B. Stability for i > 1 will also follow from 



the considerations made in section 6.2 



6 The general case 

In this section we shall prove that both solitons and black holes in which the 
gauge function w has no zeros are stable for |A| sufficiently large. Since the 
analysis required is rather involved, we shall first outline the steps in the proof, 
before describing the details. The sector with £ = has already been discussed 
in Q and section |5|, so we assume that £ > 1 in the analysis below. 



6.1 Outline of stability proof 

In the last section, we have re-formulated the perturbation equations on a 
Schwarzschild-adS and RN-adS background in terms of fully gauge-invariant 
amplitudes. Unfortunately, in the general case, this procedure is quite messy 
and we have no guarantee that the resulting equations are going to be suitable 
for analytic discussions. 

In the following, we adopt a gauge in which t = p. Somehow, this gauge 
has similar properties as r = 0, but has the advantage of being simpler compu- 
tationally. In particular, one does not lose the symmetry of the potential. For 
t = p, the momentum constraint, Eq. (15), reduces to an algebraic equation. 



V6^-^t = ^jq + ud, (24) 
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which permits one to re-express t in terms of q and d. (Note that r^p = NS is 
positive everywhere outside the horizon for black holes and positive everywhere 
for solitons.) Introducing t = p into the energy expression (12), we see that 



there are no kinetic terms involving t (or p). In particular, the kinetic energy is 
positive, and t = p is non-dynamical. The energy now has the form 

E = ^ J (^0 , tj) + {dpU , dpU) + {u , AdpU) + {u , dpAU) 

+ {U,SU)-2{U, b{d)) + n) dp, (25) 



where U = {q, g, b, c, d, e)^ and where A, S and h are obtained from A, S and h 
by removing rows and columns corresponding to t and p. The remainder terms 
TZ then are 

n = {Stt + Sii + 2Sn)t^ + 2^/l2^-i^ptq 

t b + 2V6 (u^p - u^) td + 2 {bt{a) - bi (a)) t. (26) 



Variation of these equations with respect to t, q, g, b, c, d, and e yields the same 
equations as in ( p!3|) but with p = t and where the equations for t and p are 
summed, resulting in an elliptic equation for a. 

In order to discuss the stability of the system, we proceed as follows: We 
first give initial data for t, p, q, g, b, c, d, e and their time derivatives such 
that the constraint equations t = p, ([T^), (|T7|) , ( p4|) and their time derivatives 
are satisfied. Then, we evolve the data using the hyperbolic evolution equations 
([l3|). In order to make sure that the gauge condition t = p is preserved during 
the evolution, we have to solve, at each time-step, the elliptic equation for a 
described above. In [^] we have shown that the constraint equations propagate. 
In particular, this is the case for the constraint (|2^). Next, as we have shown 
in section |3|, the terms which depend on d do not contribute to the energy E. 
As a consequence, one can check that E is conserved, i.e. dtE = 0, provided 
that homogeneous Dirichlet boundary conditions are given at the horizon and 
at infinity. Then, stability follows if we can show that E is positive. Indeed, if 
E is positive, Lij and £i cannot grow exponentially in time. In a gauge where 
the shift and the electric potential are zero, the same follows for 6gij and 6Ai 
since 



5g{t)ij = 6g{0)ij + / 2L(r)ijdr, 







SA{t)i = 6A{0)^ - J £i{T)dT. 





Since the terms involving d do not contribute, stability will hold if the oper- 
ator 

A = -dl + Adp + 9pA + 5 (27) 
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together with the remainder term TZ is positive, where 

1Z = n-2{bt{a) -bi{d))t. 

Prom now on, we therefore ignore all terms dependent on a, and shall consider 
only the energy functional 

E = E-^J[-2{U, 5(a)) + 2 (6* (a) - h (a)) t] dp, (28) 

since E is positive (and the system stable) if E is positive. 

In the first stage of the proof (subsection |6.2| ) we show that the operator A 
(which has the same form as A (27) but with slightly modified matrices A and 
S) can be written as 

A = B^B + b^b, 

where S is a first order matrix differential operator and 6 is a vector. This 
means that „4, is a positive operator. 

As a result of the factorization of A, we now have a slightly modified remain- 
der term It given by 

n = {Stt + Sn + 2Sti)t^ + 2Vl2fij^ptq 

+2V3n-futb + 2V6uJd, (29) 

(ignoring the terms containing the lapse a). Unfortunately, TZ is not positive 
(even in the Schwarzschild-adS case), so we do not immediately have that the 
energy E is positive. However, in the next part of the stability proof (subsection 



S.3), we argue that this remainder term is of subleading order compared to 
A when |A| — > oo. This is not sufficient to give positivity of E even for all 
sufficiently large |A|, since the operator A is positive and not necessarily bounded 
away from zero (this means that it may be that A could be equal to zero for 
some perturbations and then the subleading terms in TZ would be dominant, and 
may be negative). From the fact that TZ is subleading order as |A| — > oo, we can 
however deduce that E is positive for the precise value |A| = oo (in a sense to 



be made exact in section 5.3). 



It therefore remains to extend this stability for |A| = oo to sufficiently large 
values of |A|. We will do this via an analyticity argument, using the multi- 
dimensional equivalent of the nodal theorem [0] . The details of this are quite 
involved, and cover subsections 6.4-|6^. First of all, as before we ignore the 



terms depending on a in the energy functional E (p5|), and consider instead the 
modified energy functional 



E = + {dpU , d,U) + {U , AdpU) + {U , d,AU) 



+ {u, su)+n]dp. 
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We also express t (which is non-dynamical) in terms of q and d using the con- 
straint equation ([ 



where 

<Sqq = 
Sqd = 

Sdd = 

Sqb = 

<Sbd = 



fiur 

2 

.p 

^2^2^2 



V2r 



V2r 



I. The remainder term TZ then takes the form: 
Sqqq"^ + 2Sqdqd + Sddd^ + 2Sqbqb + 2Sbdbd, 



.l + 2Ar + Ar2 + ^ 



-1 + 27V + Ar^ + 



_ 1)2 



-l + 2iV + Ar^ + 



+ 



^ _^ ^P ^ """".P^T 



+ 



' ,p 
2uti r,r 



V2r„ 



71' 



2u^ 



(30) 



On varying E with respect to the perturbations 6, c, d and e, the remainder 
term will contribute additional terms to the perturbation equations. These 
additional terms are equivalent to adding the following quantities to the per- 
tinent parts of the potential S: add Sqq to ^22 in Sgrav (|8|); add Sdd to the 
((id)-entry in Syu (§); finally, add Sqd^ Sqb and Sbd to the corresponding entries 
in Sir^t (0). 

There is one additional complication. For both black hole and soliton solu- 
tions, it is found to be necessary to study the perturbation equations arising not 
from the variation of the whole of i?, but, instead we shall vary only 

E'P = E — 'P, 

where V is manifestly positive. Positivity of E is automatic if E-p is positive 
since we shall always have V > 0. Further, we can deduce that E-p > if the 
perturbation equations derived from varying Ep have no unstable modes. The 



splitting of E in this form will be seen to be necessary in section 3.5 in order 



to satisfy the technical conditions for the application of the multi-dimensional 
nodal theorem. 



In subsection we write the equations derived from the variation of Ep in 
the form 



OU 



_d_ 

dR 



+ 



-K2iR)dfU. (31) 



In equation (plD, the variable R depends on whether we are considering solitonic 



or black hole solutions and is given in subsection 6.5, and Ki, K2, L, and M 
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are real symmetric matrices depending on R. We show, in subsection that 
Ki is uniformly positive definite, K2 is positive definite, L{0) is non-negative 
and that all four matrices are smooth and uniformly bounded for all R in the 
interval [0, 00). The crux of the stability proof is to show that the operator O has 
no negative eigenvalues for all |A| sufficiently large, given that we have already 
proved that there are no negative eigenvalues if |A| = cxd. Since the matrix K2 is 
positive definite, stability follows if the operator O has no negative eigenvalues. 
Before applying the multi-dimensional nodal theorem we note that there 



are additional criteria required of the matrices L and M, which ensure that the 
essential spectrum of the operator O equals [0,oo), and that the number of 
bound states is finite. Sufficient criteria are ||T^: M — > as i? — > 00; and 
R^M + L > —(3c for some /?c < 1/4 and all sufficiently large R. We will show. 



in section |6.5| , that these additional criteria are also satisfied by our matrices L 
and M. 



The multi-dimensional nodal theorem |11] gives a method for determining 
the number of negative eigenvalues of the differential operator O, as follows. 
Fix 6 (since we have a 6-dimensional system) linearly independent real 6-vectors 
ei, . . . , eg, and ai > 0. Denote by Ua^ = [ui, . . . ,uq] the (6 x 6)-matrix whose 
columns are the solutions of the 6 initial value problems 

du ' 

Ouj = 0, ai < i? < 00, Uj{ai) = 0, Tni'^^) ~ J = 1; • • • jS. 

dR 

Then the nodal theorem reads |11]: 

Theorem: If ai > is sufficiently small and 02 > ai is sufficiently large, 
the number of zeros (counted with multiplicities) in the interval (01,02) of the 
function J := i? 1— > det Z/Z^^ (i?) equals the number of negative eigenvalues of O 
(counted with multiplicities). 

At this point it should be stressed that the function ^ must, by definition, 
vanish at oi. However, we already know that, when |A| = 00, the function has 
no zeros in the (open) interval (01,02), for sufficiently small oi. In subsection 
S.6 we show that each of the Uj (and hence 5) are analytic in A and oi in a 



neighbourhood of A = —00, and for all oi > 0. Therefore has no zeros in 
the interval (01,02) for all |A| sufficiently large (and all sufficiently small oi). 
Therefore O has no negative eigenvalues for all | A| sufficiently large, and we have 



proved stability, since the matrix K2 on the right-hand side of (31) is positive 
definite. 

6.2 Factorization of a subsystem 

Having outlined the stability proof, in this section we perform the first stage of 
this proof, by showing that the operator A is positive, where A has the same 



form as the operator A (27) but with slightly modified matrices A. and S. 

Firstly, we note that the operator A (p7|) is very similar to the corresponding 
operator in the odd-parity sector. This resemblance becomes even more striking 
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if we use the constraint equation ( p4[ ) and replace the term — 2-v/6u^ td on the 
right-hand side of equation (|^) by 

—V2fj,ju qd — 2u'^d^ . 

Redefining U = {q, —b,g, —c,d, —e)'^, the energy functional E (|2^ ) then has the 
same form as before, but where now 
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sym. 
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n = (5« + Sii + 2%)t^ + 2Vl2/i7,pig 
+2\/3^7'u th + d. 

Comparing this with the corresponding expressions in the odd-parity sector (see 
Eq. (8) of Ref. |^]) with a = 0, we see that the only difference lies in the sign 
of (6^2)23 = ('S'2)32- It is also interesting to observe that in terms of the new 
variables, the constraint equations (|l^) and (^) are 








'^dp ( -q] + ub-VX' 
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V2 



27e, 
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which, apart from the term involving t, are exactly the same as the constraint 
equations for h and b in the odd-parity sector (see Eq. (9) of Ref. [^). 

Therefore, if P denotes the projector from the even-parity amplitudes U = 
{q,b,g,c,d,e) to the odd-parity amplitudes Uodd = {h,a,b,k,c,d,e) defined in 
Ref. |], 

P -.U ^ {q,0,b,g,c,d,e), 

and if E^^ is the matrix whose entry in the column and row corresponding to c 
and d, respectively, is 1 while all other entries are 0, we have 

A = P^B^BP - 2V2fij^w {E,^ + E^,) , 

where A has the same form as the operator A ( p7| ) but with the matrices A 
and S replaced by their modified forms above. The operator B and its adjoint 
B^ were found in |^ {B is also given in Appendix]^ for completeness). Now, 
introducing also 

Q '■ Uodd = {h, a, b, k, c, d, e) ^ (0, a, 0, 0, 0, 0, 0) 
such that l-odd = PP^ + QQ^ , and noticing that 

Q^BP = 7(0, 1, 0, 0, 0, 0, 0)^ (0, 0, 0, -fi, -V2w, 0), 

we find 

A = {P^B^P) (P^BP) + 7^ (^fi'^Ess - V2fiwE^^ - "^^l^wEj^ + 2w'^E^^ 
= B^B + Fb, (32) 

where B = P^ BP and b = 7(0, 0, 0, — ^, +\/2w, 0). This shows that the oper- 
ator A is positive. Note that this factorization also works when i = 1, since in 
that case the amplitudes (g, 6, c, d) decouple from the amplitudes {g, e) so that 
the above result is still correct when g and e are removed. 

As an application of (|3^), consider the Schwarzschild-adS case, \w\ = 1: 
It can be checked from the explicit expressions in Appendix ^ that the YM 
perturbations (6, c, d, e) decouple from the gravitational perturbations. Since 
the remainder term TZ does not depend on the YM perturbations in that case, 
( |3^ shows that the YM perturbations are stable. Similarly, in the RN-adS case, 
the non-Abelian amplitudes (6, d, e) decouple, and (|32|) shows that there are no 
instabilities. In contrast, the stability of the gravitational perturbations does 
not follow from ( |3^ since in that case, the remainder term 'R, has to be taken 
into account. As a matter of fact, it turns out that TZ is not positive even in the 
Schwarzschild-adS case, and one needs to derive the master equations ( p2|) and 
( p3| ) in order to prove stability. In this sense, the gravitational perturbations are 
responsible for the stability not being topological in the even-parity sector! 

The main idea in the next section will be to show that TZ is negligible when 
|A| is large enough and to argue that in that case, the energy is still positive. 
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6.3 Stability for infinite cosmological constant 

We will now argue that the remainder term TZ (^) given in the previous sub- 
section is sub- leading order as |A| ^ oo. In order to estimate carefully the 
magnitudes of all terms in our system of equations, we firstly need to non- 
dimensionalize all quantities. All the perturbations, g, g, b, c, d and e have 
dimensions of length. Bearing in mind that Newton's constant G has dimension 
length'^, we may write, for example, 

q = qic 

and similarly for the other quantities, where £g = VG ■ Since our equations are 
linear, the factors of ic all cancel and the equations for the hatted variables are 
the same as those for the unhatted variables. Therefore we shall not distinguish 
in the following between the dimensionless or dimensionful perturbations. 
The time and radial variables are rescaled as: 

t = i(.G, P = P^G, r = flG- 

The gauge function w and all metric quantities are dimensionless. The remaining 
dimensionful quantities are: 

rh = rhiG, m{r)=m£G, u = u t^^ = 2-rf- , 7 = 7^g^- 

r iG 

We also define the positive dimensionless parameter 

L = -Ml . 

All other quantities in the pulsation equations are dimensionless. 

We are interested in the limit in which L is very large (and positive), and 
so we can introduce another positive parameter ^ = so that, equivalently, 
we may consider small ^ in a neighbourhood of ^ = 0. It is worth recalling 
at this stage that it has been proved that, in this limit, u ~ o (C~^^^) (see |Q| 
and section ^). Furthermore, we know that in this limit the metric function rh 
behaves like L, and, in turn, this means that N ^ L. We therefore define 

N = e,N = L-^N, 

and accordingly, 

7 = ?7 , 

where N and 7 will be finite as ^ — > 0. With these scalings, the potential S is of 
order as L — > 00 (since each derivative with respect to p introduces a factor 
of A^), so in order to obtain non-trivial equations we must rescale the space-time 
variables again, namely. 
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The pulsation equations now have the form 



df -dj + Adp + dpA + S -b{d) U. (33) 



where 

We also define 



AL, S = ~SL^, b = bL^. 



The pulsation equations now have a well-defined limit as ^ — > (or, equivalently, 
L — > oo). 

The remainder term TZ (where TZ is given by (p6|)) contains t terms. The 
variable t is non-dynamical and given in terms of the perturbations q and d by 



the constraint (24), which reads, in dimensionless variables: 



jq + ud] . (34) 



We therefore substitute for t to give an expression for TZ in terms of q, d and b. 
From (^^, each term in iz contains a factor of whereas the operator A is 
0(1) as ^ — > 0. Therefore, as ^ — > 0, the remainder term TZ is subleading order 
compared to the operator A. 

So far in this section we have included the terms in the perturbation equa- 



tions dependent on the lapse o. As discussed in section |6.1| , these terms do not 
contribute to the energy functional E, and so may be ignored in our study of sta- 
bility. Ignoring the b terms therefore, our conclusions apply equally well to the 



(unsealed) quantities A and TZ (see subsection 5A for the definitions of these), 
namely that, as ^ ^ 0, the remainder term TZ is subleading order compared to 
the operator A. This can also be extended to the slightly modified operator A 
constructed in subsection |6.2| , and the correspondingly modified remainder term 
TZ given in (p9|) . Therefore the (modified) remainder term TZ is subleading order 
compared to the (modified) operator A, as L — > oo. 

By the factorization in subsection we know that the operator A is pos- 
itive, and therefore we have proved that the total energy functional giving rise 
to the pulsation equations (|3^ ) is positive when ^ = (or, equivalently, L = oo), 
so that the system is stable for this value of ^. Since it is possible for A to be 
zero for some perturbations, we are not able immediately to say that the system 
is stable for sufficiently small ^ since it may be the case that for non-zero ^ the 
remainder term TZ (which is not positive) is dominant, so that the energy is no 
longer positive. The remainder of this section will be devoted to extending this 
stability to sufficiently small ^, using an analyticity argument. 

6.4 Transformation of the equations 

We now address the problem of extending our proof that the system of pulsation 
equations is stable when A = — oo to sufficiently large (and negative) A. To do 
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this, we shall use an analyticity argument based on the multi-dimensional nodal 



theorem [11] 



In order to apply the multi-dimensional nodal theorem |11|, we first need 
to cast our pulsation equations into the appropriate form. We begin with the 
pulsation equations in the form 

= {di -dl + Adp + dpA + s) u, 

where we remind the reader that A and S are anti-symmetric and symmetric 
matrices, respectively, and U = {q,g,b,c,d,e). Note that we are ignoring the 
terms in the perturbation equations which contain the lapse, since these terms 
do not contribute to the energy E. Furthermore, we derive these equations by 



firstly substituting for the variable t = p using the momentum constraint (24), 



which gives the remainder term TZ whose terms are given by (^) . Finally, 



we vary not the whole modified energy E (28), but rather E-p, which is the 



modified energy E minus a term V which is manifestly positive. This procedure 



is outlined in section 6.1, and the details for the particular cases of black hole 



and soliton equilibrium solutions are given in section 3^ below. 
Define a new perturbation vector V by 

U = RV 

where R is a transformation matrix, which may depend on p (but not t) and 
contains no derivative operators. Then the equation for V is 

= [R^Rdf - dpR^Rdp + ^Adp + dp^A + ^S] V, (35) 
where the transformed matrices ^A and are 
^ A = R^AR + ^ {R^pR - R^R,p) , 

= R^ SR -\- RF AR p — RFpAR — — (^rF R pp -\- RFppR^ . 



Since 



we set 



Aint 
int 



- T 

'A^„f 



where is a 2 x 4 matrix. If we choose F such that 

Aint — ~2"^>P' 

then ^ A = 0, and the pulsation equations simplify. With this choice, the matrix 
F takes the form 

(10 
V2 
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where the scalar function J-{p) is given by 

^(p) = f^dp = -2VG y" ^ dr. (36) 

We write the transformed potential as: 

^S = R^( 'ff ^ I'l'"* ) R, (37) 



where 



Sint SyM 



1 <? — <? —AA - ^a'^ 

^ grav — 'J grav ^-^int-^int^ 



Sint — Sint ~\~ -^int ,pi 

^Sym = Sym- (38) 

The transformation matrix R has all its eigenvalues equal to unity, but it is 
not necessarily positive definite because it may contain a Jordan block. However, 
if we introduce the vector X = (xi, X2, X3, X4, X5, xe)^, then 

X^RX = (xi + J'{p)x-if + fx6 + V2J'{p)x2 ^ 



+xl + xi + (1 - 2r\p)) xl + {l- T\p)) xl 

Therefore -R is positive definite if J^'^{p) < 1/2. Since (see section 2 and [^), 
w,r{f) ~ o(|A|~2) as |A| — > 00, for sufficiently large |A| (and an appropriate 
choice of integration constant) we can make J- as small as we like. Therefore R 
is positive definite for sufficiently large |A|. 

We now make a further transformation by introducing a diagonal matrix 
Y (depending only on p and not on t), which commutes with R. In order to 
commute with R, the elements of Y must be such that 

F = diag(yi,y2,yi,y3, 2/4,^2)- (39) 

Defining W = Y^^V , the pulsation equations ( ^5|) now take the form 

= YR^RYdfW - dp {YR^RYdpW) + "^SW, (40) 

where the new potential is 

= Y\^S -dp{R^R)Y-^Y^p-R^RY-^Y^pp\Y 

= YR^ ( 'f^r ) RY, (41) 

V "^mt Sym J 

and 



(42) 



20 _ 

^ grav — 


^ grav 


Yi ^Yi^pp, 


2o _ 
i^mt — 


iJmt 


' ^,pY2 ^Y2,p, 


'^Sym = 


^Sym 


- Y2 Y2,pp, 
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and we have defined diagonal matrices Yi and Y2, which form the blocks of the 
matrix Y, as follows: 

Yi = diag(yi,y2), 
Y2 = diag(yi,y3,y4,y2)- 
At the moment we leave the functions in the matrix Y to be quite general. 



They will be fixed in the next subsection p.5| , when we discuss the form of the 
equations in detail for black holes and solitons in turn. 

For the moment, we write the pulsation equations ( |40| ) as 

- K{p)d!W = -dp {K{p)dpW) + ^SW, (43) 

where 

K = YR^RY. 



Equation (43) should be compared with the form (|3l| ) required for the application 
of the nodal theorem The first issue is whether our coordinate p has the 
required domain [0,oo). This will be addressed in the following subsection for 
the black hole and soliton solutions. For the time being, we assume that we can, 
if necessary, change to a coordinate R (which will be a function of p) which does 
have the required range of values. If 

then the equations (^) become 

which is of the form (^), with Ki = KQ and K2 = KQ~^. In our work in the 
following section, the function Q will always be of one sign, and so we consider 
instead (by multiplying throughout by -1 if necessary) 

- K\Q\-'d^W = - A (kIQI^) + '~S\Q\-'W, (45) 

so that the matrices Ki = K\Q\ and K2 = K\Q\^^ are positive definite. As dis- 
cussed in section |6.1| , in order to prove stability we require that Ki is uniformly 
bounded and uniformly positive definite, a condition which will be considered in 
the next subsection. In order to apply the generalized multi-dimensional nodal 
theorem, it remains to show that the matrix ^S^|(5|~^ can be written in the form 

^S\Qr' = ^L{R) + M{R), (46) 

where L and M are uniformly bounded matrices on i? G [0, 00), and L{0) is non- 
negative. Furthermore, we also require M ^ as i? ^ cx), and MR^ + L > -f3c 
for some f3c < 1/4 and all sufficiently large R. This is the subject of the next 
subsection. 



30 



6.5 The nodal theorem matrices 

Having cast our pulsation equations into a form (|4^) which wiU ahow us to apply 
the (extended) multi-dimensional nodal theorem, we now have to verify that the 
matrices in (p5| ) do in fact satisfy the conditions of the nodal theorem. Namely, 
-f^lQI and -f^jQI"^ must be uniformly bounded and positive definite, and the 
potential '^S\Q\^^ must have the form (^). In addition, M ^ as i? ^ oo, 
and MR^ + L>-(3c for some /?c < 1/4 and all sufficiently large R, which will 
be the case if ^S\Q\-^ -^OasR^oo, and ^S\Q\-^R^ > -/3c for all sufficiently 
large R. So far we have also left open the precise identity of the coordinate R 
for use in the application of the nodal theorem. 

In this subsection we shall address these issues. The two cases, black hole 
and soliton solutions, need to be addressed separately because of the different 
boundary conditions on the equilibrium solutions. The black hole perturbation 
equations are more readily cast into an appropriate form, so we deal with them 
first. 



6.5.1 Black holes 

For black hole equilibrium solutions, the usual "tortoise" coordinate p has the 
domain {—oo, pmax], tending to — oo close to the event horizon, and approaching 
the finite value pmax at infinity. By a suitable choice of integration constant in 
the definition of p, we may set pmax to be equal to zero. Therefore p S (— oo, 0]. 
Our application of the multi-dimensional nodal theorem requires a variable R 
having domain [0, oo), so we define 

R = -p. 

It should be noted that R = corresponds now to infinity, and i? — > oo corre- 
sponds to the event horizon. 

In this situation, the quantity Q (see (|4^)) is equal to —1, and so \Q\ = 1. 
We set the transformation matrix Y (p9|) to be equal to the identity, so that the 



perturbation equations are now of the form (45) 



where 



K = R^R, 



and is given by (|37|). 

The first necessary condition is that the matrix K is uniformly bounded 
and uniformly positive definite on G [0,oo). Recall that R depends on the 
function J-{R), defined by (|3^). We are interested in those black holes for which 
the gauge function w has no zeros, in which case, from the field equations, the 
function w is a monotonic function and so w^r is of one sign. To see this, suppose 
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that w > 0. First note that the field equations (0) imply that w cannot have a 
minimum in the interval < u; < 1 nor a maximum if w > 1 [12|. Furthermore, 
a point of inflection is possible only if t;; = 1, and it was proven in |^] that 
for solutions in which w has no nodes, the function w cannot cross ±1. Taken 
together, these statements mean that w can have no stationary points, and this 
argument applies equally well to soliton solutions or if w < everywhere. 

Given the asymptotic behaviour of w at inflnity (see section |2|) , in this case 
the integral defining JF (^6|) is convergent at infinity (which corresponds to -R = 
0). This means that the matrix R is uniformly bounded on i? G [0, cxd), so 
that K is also uniformly bounded. It is also straightforward to show that K is 
uniformly positive definite on R £ [0, oo) in this case. 

Since the matrix R is uniformly bounded and positive definite, to show that 
the potential can be written in the form (^6|), it is sufficient to show that this 
is the case for the potential matrix 

1 Q 1 O 

grav >-'int 



whose components are given in (pq). Remember that R = (which is where the 
potential can be singular) corresponds to infinity. It is simplest to consider each 
part of the potential separately. 

Firstly, from (|3^), there are terms in the potential proportional to Aj„f 
and/or its derivative with respect to p. These depend only on the gauge func- 
tion u ( |To|) and its derivative, respectively. Both these functions are regular 
everywhere outside the event horizon and uniformly bounded for R E [0, oo). 
Therefore those parts of the potential containing Aint or its derivative may be 
included in the matrix M (46). 

It now remains to consider only the terms Sgrav , Sint and Sym , remembering 
to include the additional contributions from the remainder term (^0|). Careful 
analysis of each component in these matrices reveals that all components are 
regular and uniformly bounded everywhere outside the event horizon (so that 
they may be included in the matrix M) with the sole exception of the quantity 

r 

which arises in the ^33 component of the Sgrav matrix (see (|8|)). At infinity 
(-R — > 0), we have 

We therefore define a matrix L which has as its only non-zero element 

J- _ ' ,pp 

r 

Near the event horizon i? ~ — log(r— r/^) — > cxd, and so L33 0. Thus the matrix 
L is symmetric, regular and uniformly bounded on i? € [0, cxd). In addition. 
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-^33(0) = 2, so that the matrix L{0) is non-negative, as required. We comment 
here that a crude estimate of the behaviour of 

„2 



r 



a 



which arises in the component of Sgrav suggests that this is divergent at 
infinity. However, a precise calculation, using the asymptotic forms of the metric 
functions (section |2|) shows that in fact this function is regular everywhere outside 
the event horizon. 



Having written the potential in the required form (46), it remains to show 
that this potential vanishes sufficiently quickly as i? ^ 00, that is, 



for all sufficiently large R, where /3c is a constant which is less than 1/4. Recall 
that R ^ 00 corresponds to the event horizon of the black hole. As the event 
horizon is approached, all terms in the potential tend to zero at least as 
rapidly as (r — r/i)^/^, with the exception of the following terms: 



-'qq 



2rh 



+ 0{r - Th) 



S- 



22 



-S\rh)Nl.{rh) + 0{r 



where Sqq is the extra contribution to 5*22 in S 



grav 



given in (30). In addition. 



the (66)-component of Sym (^) has the same behaviour as 822- Therefore there 
are two entries in the potential matrix which approach a positive constant as 
R ^ CO (since N^ri^h) > 0). We deal with these by defining a manifestly positive 
quantity V: 



V 



2rh 4 ' 



1 



which we subtract from the modified energy functional E before varying Ej> — 
E — V to yield the perturbation equations. These equations are identical to 
the previous ones apart from two terms in the untransformed potential S: the 
(g(?)-term in S gfav ^ which, has the positive constant 

2rh 4 

subtracted, and the (bb)-term in Sym, which has the positive constant 

subtracted. This does not affect any of the analysis so far in this subsection, in 
particular the definition of the matrix L is unchanged, and the matrix M simply 
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has appropriate positive constants subtracted from its (qq) and (bb) entries. 
Therefore both L and M remain uniformly bounded for all R G [0, oo). 

However, now all entries in the potential vanish at least as quickly as 
(r — r/i)^/^ as r ^ r/i (or, equivalently, R oo). Since, near the event horizon, 
R ~ — log(r — r/i), this means that the potential is vanishing like e~^/^ as 
i? — > cxD. Therefore ^SR^ — > as ^ cxd, so it is certainly the case that 

^SR^ > -Pc 

where Pc < 1/4 for all sufficiently large R. 

Therefore, for black hole equilibrium solutions, we have constructed matrices 
satisfying the conditions necessary for the application of the multi-dimensional 
nodal theorem. 



6.5.2 Solitons 

The corresponding analysis for soliton solutions is rather more complex. We shall 
begin by defining V = {'y'^fi^r'^/2r'^p)q'^ as it is manifestly positive, and consider 
the perturbation equations derived from the variation of E-p = E—V (as outlined 
in section |6.lD . This makes some of the algebra later in this subsection more 
tractable. We can prove stability for soliton solutions if these pulsation equations 
have no unstable modes. Therefore, for the pulsation equations for solitons, we 
remove a term 

from 

in S grav All other terms in the potential are unchanged. 
For soliton solutions, the "tortoise" co-ordinate p has a finite range of values, 
so we shall instead begin with the co-ordinate r, which has the domain [0, oo). 
We define a new co-ordinate R by 

R = r~^, (47) 

which also has the domain [0, oo). However, R = now corresponds to infinity 
(r — > oo), while R ^ oo corresponds to the origin (r = 0). This mapping of 
infinity to i? = can be motivated from our experience in dealing with the 
black hole solutions in the previous subsection, as we already know in detail the 
behaviour of the potential at i? = 0, which is crucial for applying the multi- 
dimensional nodal theorem. The specific choice of co-ordinate ( ^7| ) will be seen 
below to enable us to put the potential in such a form that the precise conditions 
necessary for the application of the nodal theorem are satisfied. 
With R defined by (|^, the quantity Q (||) is given by 

Q = Ar"^NS. 

Since this is negative, we shall use |Q| in the perturbation equations (as in (^5|)). 
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As with the black hole solutions in the previous subsection, T ( pq ) is regular 
and uniformly bounded on the whole of r G [0, oo) and therefore so too is the 
transformation matrix R. 

For our stability proof, we require the matrix -K^IQI to be uniformly bounded 
and uniformly positive definite. However, \Q\ ~ r^l"^ as r — > oo, while R is 
uniformly bounded, so we need to introduce a second transformation matrix Y 
(see section such that -f^|Q| = YliF T{Y\Q\ is uniformly positive definite 
and bounded. Let 

which satisfies the requirement that Y should commute with R. The factor 
1/2^/2/i will be seen later in this subsection to ensure that the behaviour of the 
potential as i? ^ oo satisfies the conditions of the nodal theorem. Then 

KU = ^fl-fl 

which is uniformly bounded and positive definite on i? G [0, oo), as required. 

With the additional transformation given by the matrix Y ^ the potential 
is now ^SIQI"^, where is given by (|4l|). Since the matrix R is uniformly 
bounded and positive definite, as in the previous subsection, we only need to 
show that the matrix 



\Q\-'Y 



2 a 2 Q 

l-^ I grav i^int 



can be written in the form (^). The elements of this matrix are given by (|42|). 
As with the black hole solutions, we consider the various contributions to the 
potential in turn. 
Firstly, 

as r ^ {R ^ oo) and 



Y = , r\T + 0(rx) 



Y = r-\l + Oir-h 

as r — > oo (i? ^ 0), where we have used the behaviour of the metric functions 
in the asymptotic regions (see section ^). Therefore, as r ^ 0, 

Y{Y-^Y,p)Y\Q\~^ = -^r^SiOy'l + Oir^, 

Y{Y-'Y,^^)Y\Q\-' = --l^rl + 0{r'y, (48) 
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and as r — > oo, 

Y{Y^^Y,,,)Y\Q\-^ = -^rl + Oir-^). 

Using the fact that the matrix F and its derivative with respect to p are both reg- 
ular and uniformly bounded, for R,r G [0, oo), the contribution to the potential 
F^pY2^Y2 ,p is also regular and uniformly bounded, and so can be included in the 
matrix M. However, as r ^ oo (i.e. R 0), the matrix Y {Y'^Y^pp) Y\Q\-^ 
diverges like r ~ R~^. Therefore this part of the potential will need to be 
included in the L matrix, which we will consider in detail below. 

Next we consider the contribution due to the R transformation, given in 
(|38|). The matrix Aj„( and its derivative Aint,p are both regular and uniformly 



bounded on r E [0, oo), and so those terms in the potential (38) depending on 
Aint can be included in the matrix M. 

It remains to consider the following part of the potential 



\Q\-'YSY = \Q\-'Y( -^T ^ 
whose components are given in section ^. From our analysis of the potential for 



black hole solutions (subsection |6.5.1 ), we know that the potential component 

^grav 



S33 in Sgrav diverges like as r — > 00, all other terms in S being regular. When 



multiplied hyY^\Q\-^, this gives a divergent term of the form 
^. + 0(.-') = ^fl- + 0(/i^). 



Therefore we define a matrix L (|4q ) to be equal to + '^L, where the only 
non-zero entry in is 

1 



1 ^ 

2^2' 



and 



L33 

2 



Therefore the matrix L is uniformly bounded on i? G [0, 00), and, furthermore, is 
constant. Since fj? > 0, we also have that L is positive definite, so in particular 
L{R = 0) is positive definite, as required by the nodal theorem. 

We subtract LR~^ from the full, transformed potential ^^IQI^^ to give the 
matrix M, which is then uniformly bounded for R £ [0, 00). The final criterion 
we need to check before the multi-dimensional nodal theorem can be applied in 
the next subsection is that the potential 25^|Q|~i vanishes sufficiently quickly as 
i? — > 00 (i.e. r — > 0) that 

'S\Q\-'R' > 
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for all sufficiently large R, where Pc is a constant less than 1/4. This requires a 
detailed study of the behaviour of the potential near the origin r = 0. 

We have already considered the form at the origin of the contributions to 
the potential due to the Y transformation (^). Those parts of the potential 
( |38|) arising from the R transformation, as already observed, are bounded at the 
origin, and so when multiplied by ~ as r ^ 0, they vanish at least 

as quickly as R~^, and so need not be considered further here. 

The behaviour of the untransformed potential S at the origin r = is com- 
plicated due to the presence of 7 and its derivatives. Near the origin, 



7 



5(0) 



0(1). 



By considering each term in the potential in turn, we arrive at the following 
explicit expression for the untransformed potential near the origin: 



S\0) 



( 2 



V 

M+0{r 



- A 


-2^/\ 











2y/\ 


A 

















A + 4 


2 


-72/2 








2 


A + 4 














A + 4 








-\/2A 


-^/2A 






2A 
2A 




A 



+ 0(r-i) 



where, in the second line above, we have defined the matrix M. Notice that the 
entry in the first row and column in the matrix M which is negative for £ > 2 
originates from the (g(;)-term in the potential ^S, including contributions from 



both the untransformed potential term 5*22 (§) and the remainder term Sqq (|3_ 
Multiplying this by and adding the contribution to the potential from 

the Y transformation (whose behaviour near the origin is given by (^)) 
have, for the complete potential '^S, 



we 



1 3 

4/^2 64^2 



-j^M + — 

4/x2 64yu2 



I]R-^ + 0{R- 



as R ^ 00. Therefore the potential vanishes like R^^ as R ^ 00. 

In order to check whether '^S\Q\^^ R'^ > —j3c for all sufficiently large i?, where 
f3c < 1/4, we need to calculate the eigenvalues of the matrix Ai. Firstly, arising 
from the bottom right hand part of the matrix, 

£{£-!), {i + l){i + 2), 

both these eigenvalues repeated twice. For i > 1, these eigenvalues are non- 
negative. From the top left-hand part of the matrix we have the single eigenval- 
ues 

-A, \ + 2 = ii^. 
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One of these eigenvalues is positive, the other negative, leading to a negative 
eigenvalue oi'^S\Q\~^R'^. However, this negative eigenvalue of'^S\Q\~^R'^ takes 
the form 

_1 _35___ _1 

Therefore '^S\Q\~^R'^ > —(3c for all sufficiently large R, where we have defined 
the constant /3c < 1/4 in the equation above. 

As an aside, we note that the eigenvector corresponding to the negative 
eigenvalue of M has the form 

(g = ^/A5,5,0,0,0,0). (49) 

We are working in a gauge in which t = p. This condition does not completely 
specify the gauge, and some freedom remains in our quantities q,g,b,c,d,e. 
Physically we are interested only in gauge-invariant quantities. One such gauge- 
invariant is (, given in (^9|). Since we are working in a gauge in which t = p, 
near the origin has the form 

Therefore the eigenvector ( |4^ ) corresponds to a perturbation in which the gauge- 
invariant C vanishes near the origin, and so represents a pure gauge perturbation. 
The reader may wonder, at this point, whether the existence of the negative 
eigenvalue in A4 is due to our leaving part of the potential term 5*22 in the 
remainder V. In fact, even if we had not performed that step, we would still 
obtain a negative eigenvalue in the corresponding Ai matrix. However, in that 
case the eigenvalue is considerably more algebraically complicated than the value 
—A we have in this case, and, furthermore, the associated eigenvector is not so 
easily interpreted as being pure gauge. 

We note furthermore that, in order to construct the potential in such a way 
that the conditions for the application of the nodal theorem are satisfied, it was 
essential that we used a co-ordinate R which was zero at infinity (r — > oo), and 
tended to infinity at the origin (r = 0). Otherwise, we would have been forced to 
construct an L matrix such that, at i? = 0, this matrix had a negative eigenvalue, 
whereas the nodal theorem requires that L is positive definite at -R = 0. By 
setting i? — > cxD at the origin r = 0, the matrix having the negative eigenvalue is 
instead the limit of as R ^ oo, which need not be positive definite, 

but instead must simply be bounded below by a constant strictly greater than 

— 1/4. Including the factor of (2\/2/u)~^ in Y, we have been able to ensure that 
the only negative eigenvalue of this matrix at ^ (X) is strictly greater than 

— 1/4, thus satisfying the criteria. 

To summarize, in this subsection we have shown that, for soliton solutions, 
the potential has the form which enables the application of the multi-dimensional 
nodal theorem. We are now in a position to complete our stability proof. 
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6.6 Analyticity argument 

In the previous subsection we cast the pulsation equations in a form which 



enables us to use the multi-dimensional nodal theorem. From section 6.3, we 
know that there are no unstable modes of the system when |A| = oo. We now 
wish to extend this result to sufficiently large |A|. To do this, we consider the 



function 5^ := i— > detUai{R) (see section 6.1 for definitions), which we know 
has no nodes in the interval (01,02) for sufficiently small ai and sufficiently 
large 02 when |A| = 00. In this subsection we shall show that ^ is analytic in 
|A| (or, equivalently, ^ = — A~^), and therefore, in a neighbourhood of ^ = 0, 
the function ^ will also have no zeros (in an appropriate interval (ai, 02))) which 
means that our system of pulsation equations has no instabilities. 

Since the function ^ is defined in terms of the matrix Uai , having columns 
Uj{R), where the Uj are the solutions of the initial value problems: 

duj ' 

Ouj = 0, ai < R < 00, Uj{ai) = 0, "777(^1) = ^j) J = 1; • • • )6 

dR 

(50) 



(where the operator O is defined in (|3l|)), ^ will be analytic if we can prove that 
each Uj{R) is analytic in ^. Prom the background solutions are analytic in ^ 
in a neighbourhood of ^ = 0, so all the elements in the matrices in O are analytic 
in ^ for sufficiently small ^. The differential equation for Uj has regular singular 
points at i? = and 00, but all i? G [ai, 02] with ai > and 02 > fli are regular 
points of the differential equation (50). 

Define a new independent variable Rhy R = R — ai, and make R and ^ into 
dependent variables. Then we consider the differential equations in (|50| ) (but 
now with R as the variable) together with 

^ = 1 ^ = 
dR ' dR ' 

The initial conditions are now: 

du ' 

Uj{0) = 0, -4(0) =e,-, R{0) = ai, m = ^■ 
dR 

With this set-up, i? = is a regular point of the system of differential equations, 
and so standard theorems (see, for example, proposition 1 of |jl2|) tell us that 
we have a solution of the initial value problem, at least in a neighbourhood of 
R = 0, which is analytic in R and the initial parameters oi and ^, for all oi > 
and ^ sufficiently small. Furthermore, since the differential equations are linear, 
and the only regular singular points are at i? = — oi and R = 00, any solution 
which exists in a neighbourhood of i? = can be extended to R sufficiently large 
and positive, and will remain analytic in R, ai and ^. 

To summarize, for all ai > and sufficiently small ^, we have solutions Uj to 
the initial value problems ( |50| ) which are defined on S [oi, 02] for all 02 > aij 
and are analytic in ^ and R. Therefore, for all oi > and sufficiently small ^ we 
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have a function 5 which is defined on i? G [01,02] for all 02 > ai and analytic 
in ^ and R. We have already shown that, when ^ = 0, the function ^ has no 
zeros in the interval (01,02) for all sufficiently small oi and sufficiently large 
02. We stress once again that 5 rnust vanish at = oi (i.e. = 0) by the 
initial conditions on the Uj (^0|). In other words, when considered as a function 
of R, the function ^ has no zeros on G (0, 00), for sufficiently small oi, when 
^ = 0. We now invoke analyticity to state that for sufficiently small ^ and oi, the 
function ^ will still have no zeros on R £ (0, 00), i.e. in the interval R G (oi, 02) 
for sufficiently large 02. Therefore the system of pulsation equations remains 
stable for sufficiently small (non-zero) ^. 



7 Conclusions 

This paper completes the proof that solitons and black hole solutions to 5u(2) 
EYM theory with a negative cosmological constant are stable if the magnitude 
of the cosmological constant is sufficiently large. Here we considered the even- 
parity sector of perturbations, which proved to be considerably more complex 
than the odd-parity sector studied in |Q. It is testament to the power of the 
curvature-based formalism employed here that we have been able to complete 
this proof. Unlike a metric-based formulation, this formalism yields a symmetric 
wave operator for the perturbations. Using the nodal theorem, the idea is then 
to count the number of nodes of the 'determinant of the zero modes' which gives 
the number of bound states of the spatial part of the wave operator. Here, we 
have shown analytically that this determinant has no nodes by taking the limit 
A — > — 00. Using the fact that the background solutions depend analytically on 
A we have been able to extend this result to large |A|. However, by counting 
the nodes using a numerical code, it would be straightforward to extend the 
stability analysis to any value of A. Furthermore, using a numerical code, we 
also expect that the non-spherical stability of the asymptotically flat black holes 
and solitons |§, |3| and of the topological black holes [13| could be discussed on 



the same lines as in the present article. On the other hand, our analysis also 
revealed the limitations of the curvature-based formalism: Since one does not fix 
the residual gauge completely, there can still be gauge modes. For solitons, we 
have seen that these modes make it difficult to write the perturbation equations 
in the form required by the nodal theorem. 

Physically, our result shows that the presence of a large, negative cosmo- 
logical constant stabilizes black hole hair, and also non-Abelian solitons. This 
is important in the light of the recent revival in the study of asymptotically 
anti-de Sitter spacetimes, due in part to the (rather less recent) idea of "holog- 
raphy". The concept of holography (and its more recent relative, the adS/CFT 



(conformal field theory) correspondence, see, for example, |19] for a review) cen- 
tres around the idea that behaviour (for example, of a supergravity theory) in 
the bulk of an asymptotically adS spacetime is governed by the behaviour of a 
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related theory on the boundary. In terms of these ideas, the black holes con- 
sidered in this paper are manifestly "non-holographic", since their structure is 
not determined by quantities (such as magnetic charge) measurable at infinity. 
The consequences of the existence of these stable black holes for the adS / CFT 
correspondence remain to be investigated. 

On a classical and semi-classical level, there are various open avenues of 
study. Firstly, we would like to extend our current work to seek rotating black 
holes and solitons in this model. If the rotation is slow, one could use the lin- 
earized theory to find such solutions (they are linked to the zero modes of the 
system considered here). In the asymptotically flat case, slowly rotating solu- 



tions were found in ||2^ and some of their rapidly rotating counterparts in |21]. 
Secondly, the behaviour of quantum fields on hairy black hole backgrounds has 
received little attention to date (no doubt partly due to the classical instability 
of most hairy black holes). These questions shall be the subject of future work. 
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A On the stability of the RN-adS solution 

In this appendix, we discuss the positiveness of the operator 

-3Af ^/4GA 



A+ = -dl + iV 



\/4GA 3M 



where the functions A^, \J and W are given in section p^ . 

We start by showing that the function /r(?^) = A/2 + ?>M/r — 2G/r^ which 
appears in the denominators of U and W is positive for all A > and all 
f ^ fh- Introducing the dimensionless quantities x = r/M, = G/M"^ and 
I = -3/(AM2), we can write as 



iV(^) = l-- + ^ + T 

X X'^ I 



Using the fact that N^x{xh) > and N{xh) = 0, one can show that 3xh — > 
x| > 0. Therefore, the function 3/x — Iq^ jx^ is positive outside the horizon, 
and friy^ > for all r > rh ■ 

When i > 2, one can generalize Chandrasekhar's arguments [22| in order 
to relate the operator A-\- to the corresponding operator A- in the odd-parity 
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sector by a supersymmetric-like transformation. This works as follows: First, 
we diagonalize the operator A-\- , obtaining the two decoupled operators 

A+ = -d^ + N[U + aiW] 

where ai = -VOM^ + AGX , a2 = +\/9M2 + 4GA . One can check that A+ can 
be factorized according to 

Ai+ = BiBl-Lol 



where 

^i = Op + -7T — ■ r + , B\= -dp + — — ■ + uji , 

r[Xr + Qi) r{Xr + qi) 

with qi = 3M — ai and where uji = {£ — l)i{£ + l){£ + 2)/2qi are the algebraic 
special frequencies. (Note that (Ar + qi){Xr + 52) = ^Xr'^fr > 0, so Ar + 52 is 
positive.) The supersymmetric partners of A+ are 

Explicitly, one finds 

where C/_ = {£{i + 1) - 3M/r + AG/r'^)/r'^ and = which is equivalent 
to the operator in the odd-parity sector! Since we have shown in |^] that there 
are no unstable modes in the odd-parity sector, and since for perturbations with 
compact support, A+ and A- have the same spectrum, it follows that there are 
no unstable modes in the even-parity sector either. At this point, it is interesting 
to note that the functions 



e 



Ui 



A + f 

fulfill B\ui = and therefore are eigenfunctions of the operators AiJ^ with neg- 
ative eigenvalue —iof. However, these functions do not satisfy homogeneous 
Dirichlet boundary conditions at infinity and should, therefore, not be consid- 
ered in our stability analysis. Nevertheless, the fact that Ui are finite as r ^ 00 
could affect the equivalence of quasi-normal modes in the odd- and even-parity 
sector. It has been shown recently |17] that - unlike for the Schwarzschild or 
RN case - the quasi-normal frequencies of the Schwarzschild-adS solutions are 
different in the odd- and even-parity sector. 

Finally, we discuss the case £ = I, where A+ reduces to 

-d^ + N[U + 3MW] . 

Explicitly, one finds 

(2/^)2 [U + 3MW] = [9lx^ + Aq^x^ - ISlq^x'^ + 16//x - 4/^^] . 

Using the fact that q^ <1 for a black hole, it is not difficult to show that this is 
positive when x > 1. This shows the stability of all black holes with > M. 
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B Factorization of the spatial operator in the odd- 
parity case 

In the odd-parity sector ||8|, we were able to factorize the spatial operator such 
that the perturbation equations assume the form 

Uodd + B^BUodd = 0, 

where Uodd is the vector containing the odd-parity perturbations. For £ > 2, the 
operator B is given by 



B = d,+ 



Ci cf 

Ct D{Xo + T{p)X,)D-A2 



where 



V 



2 I L 









u 


\ 


7 









-%£. 

7 


/ 



/ -Vx 





V 







— V 

w 





JL. 

V2 



In addition, D = diag(/i\/A, /u, V^, /iy^ A/2) and the matrices Xq, Xi and 
A2 are 



Xr 



Xi 



/ —Xn'^^-^ + fiwuv sym. 

—fiw u 2ww^p 

fiu -2'W,p 

\ 210"^ (1 — w'^)u 2w'^w^p 

I 2uP'v^ sym. sym. sym. \ 
—2nP'v 2w'^ sym. sym. 
2wv -2w 2 

V f2WV -f2W f2 

/ 1 \ 



sym. 
■sym. 


21010 , 



sym. \ 
sym. 
sym. 

-f2WW^p j 



sym. 
hfl J 



u 
72 






V -1 / 



where fi = X + 2w'^ and /2 = — 2w'^. 

Finally, the function T{p) has to satisfy the differential equation 



- 5„r + AT^ +BT + C = 0, 



(51) 
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with 



^^{X + 2w^) + l 



'^(^X + 2w^)^ + 2X + W 



2 2 \ / , 2 \ 2 

- M A ( - + U- ) 7 



In we have shown that equation (51) admits global solutions, at least when 
|A| is sufficiently large. (There are two small errors in the expressions for B and 
C which are published in Ref. [^]. However, these errors do not affect our proof 
that there exist global solutions to equation (^I[).) Following the same lines, it 
is also easy to show that equation (^) admits global solutions when \w\ = 1 or 
w = 0. 

In the sector with i = 1, the matrix- valued operator B can even be given 
explicitly, see [Q. 
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